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Abstract. We calculate for La 2 Cu04 the phonon- induced redistribution of the 
electronic charge density in the insulating, the underdoped pseudogap and the more 
conventional metallic state as obtained for optimal and overdoping, respectively. The 
investigation is performed for the anomalous high-frequency-oxygen-bond stretching 
modes (OBSM) which experimentally are known to display a strong softening of the 
frequencies upon doping in the cuprates. This most likely generic anomalous behaviour 
of the OBSM has been shown to be due to a strong nonlocal electron-phonon interaction 
(EPI) mediated by charge fluctuations on the ions. The modeling of the competing 
electronic states of the cuprates is achieved in terms of consecutive orbital selective 
incompressibility-compressibility transitions for the charge response. We demonstrate 
that the softening of the OBSM in these states is due to nonlocally induced dynamic 
charge inhomogenities in form of charge stripes along the CuO bonds with different 
orbital character. Thus, a multi-orbital approach is essential for the CuO plane. The 
dynamic charge inhomogeneitics may in turn be considered as precursors of static 
charge stripe order as recently observed in La2- £I ;Ba :r Cu04 in a broad range of doping 
around x = 1/8. The latter may trigger a reconstruction of the Fermi surface into 
small pockets with reduced doping. We argue that the incompressibility of the Cu3d 
orbital and simultaneously the compressibility of the 02p orbital in the pseudogap 
state seems to be required to nucleate dynamic stripes. 
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In recent work [TJ |2] we have developed a description of the electronic charge 
response and phonon dynamics for p-doped and n-doped cuprate high-temperature 
superconductors (HTSC's). This approach is based upon a microscopic modeling 
of the change of the charge response across the insulator-metal transition in terms 
of the compressibility of the electronic system as a primary tool to characterize the 
corresponding competing ground states depending on doping. This is achieved starting 
from the insulating state via a metallic pseudogap state to a more conventional metallic 
state by consecutive orbital selective incompressibility-compressibility transitions in 
terms of strict sum rules of the charge response given in Ref. [1] and listed below. 
These sum rules are a proper instrument to describe the competition between electronic 
localisation and itinerancy which is at the heart of the physics in the cuprates. 

In the conventional metallic state the electronic partial density of states (PDOS) 
at the Fermi level Z K (ef) is related to the polarizability matrix in the longwavelength 
limit (q — > 0) according to 

^II KK ,(q 0) = Z k (e f ) , (1) 

and the total density of single particle states at energy e is given by 

Z{e) = Z « ( £ ) • ( 2 ) 

K 

the orbital degrees of freedom in the elementary cell of the crystal. In our 
calculation for LaCuO, Cu3d, Cu4s and 02p orbitals in the CuO plane are taken into 
account, q is a wavevector from the first Brillouin zone (BZ). On the other hand, for 
the insulating state we obtain the sum rules 

^n KK ,(q 0) = O(q) (3) 
k' 

and 

(q 0) = O (q 2 ) . (4) 

The sum n KK / (q — > 0) is equal to p 2 K with p the average density and K the 

compressibility of the electronic system. The latter provides a measure of the charge gap 
in the electronic spectrum because K vanishes as a function of the chemical potential 
in the gap region. Moreover, 1/K provides a measure of the change of the chemical 
potential with particle number. 

Approaching the delocalization-localization transition from the conventional 
metallic region representative for optimal and overdoping, when p-doping is decreased 
the Cu3g? component of the wave function is admitted to become incompressible, 
insulator-like in the underdoped state according to our modeling. This means that 
the correlated 3d electrons are allowed to become more localized as compared with the 
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conventional metallic state and disappear from the Fermi surface (FS). In terms of the 
sum rule we then have 



Accordingly, the partial density of states is suppressed at the Fermi level for the 
correlated, localized Cu3<i orbitals (orbital intrinsic incompressibility due to correlation) 
but not for the more delocalized 02p orbitals at the O xy sublattices in the CuO plane, 
where the holes are predominantly injected in the p-type cuprates. So, we have an 
orbital selective compressible, metallic charge response only for the 02p states with a 
renormalized PDOS Z K (e f ) and a loss in the total density of state at Bp (pseudogap) 
because of the incompressibility of the Cu3d states according to equation The loss 
in spectral weight also is never recovered on entering the superconducting state. The 
incompressible regions in the pseudogap state compete with overall metallic behaviour 
and with superconductivity 

Our modeling of the pseudogap state in terms of compressible, metallic 02p states 
and Cu3g? states treated as incompressible and disappearing from the FS because of 
strong correlation is supported by calculations performed within the self-interaction- 
corrected (SIC) local spin density (LSD) approximation for La 2 Cu0 4 [3]. Here it is 
shown that the self-interaction pulls the occupied localized Cn3d states below the 02p 
band, only little Cu3<i admixture is found in the 02p band. 

Furthermore, our modeling of the pseudogap state corresponds with a smaller Fermi 
volume as compared with the conventional metallic Fermi-liquid state because the Fermi 
volume is determined by the density of the compressible metallic charge carriers (holes) 
of dominantly 02p-type alone while the Cu3d states are incompressible as in the Mott 
insulator and cannot contribute to the Fermi volume quite in contrast to the conventional 
metallic state where both the 02p and the Cu3<i states are compressible metallic leading 
to a large Fermi volume. The pseudogap state is separated from the conventional overall 
compressible Fermi-liquid state by an orbital selective compressibility-incompressibility 
transition of the Cu3d orbital. In the context of our modeling of the pseudogap state 
it is interesting to point to another modeling of this state within the framework of 
the so called fractionalized Fermi liquid. For a recent discussion see e.g. Ref. [I] and 
references therein. In this approach an idea is important which is based on the concept 
of an orbital-selective Mott transition. For the latter, e.g. in a two band model, we 
have one metallic band with Fermi-liquid like quasiparticles while the other band, e.g. 
related to localized d electrons, undergoes a Mott transition to an insulating state due 
to strong electronic correlations. Such a transition is analogous to the compressibility- 
incompressibilty transition of the Cu3d orbital in our approach. 

Because particle density fluctuations scale with the compressibility of the system 
the sum rules introduce varying fluctuations of the particle number in the different 
metallic states of the HTSC according to our modeling. The tuning parameter for these 
fluctuations of the particle number is doping and the pseudogap state may be considered 




Cu3d 
02p . 



(5) 
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as a state with reduced particle density fluctuations as compared with the optimal and 
overdoped state. Moreover, the partial incompressibility of the Cu3d orbitals competes 
with superconductivity because the phase fluctuations A0 of the order parameter grow 
with decreasing particle fluctuations AN due to the uncertainty relation AN A0 ~ 1. 

Qualitatively, the incompressibility of the Cu3<i states also corresponds to a 
compartmentalisation of configuration space. This means that some parts of direct 
space cannot be approached or are hardly accessible to the charge carriers. On the other 
hand, this should be accompanied by a reciprocal compartmentalisation of momentum 
space, according to the uncertainty relationship Ax Ak ~ 1. As a consequence the 
momentum should have corresponding fluctuations and the FS of the conventional 
metallic state becomes fuzzy when passing to the underdoped state and possibly 
reconstructs. ARPES studies and quantum-oscillation measurements seem to reflect 
this in form of a reconstruction of the large FS in the optimal and overdoped state into 
small Fermi surface pockets with Fermi-liquid like charge carriers in the underdoped 
pseudogap state characterized by a loss of density of states at the Fermi level consistent 
with our modeling. 

In our orbital based local picture the pseudogap state of the cuprates looks like a 
"two component" electronic structure where a real space organisation of the low lying 
charge excitations is achieved in form of a metallic charge response by mobile holes on 
the oxygen network in the CuO plane that is blocked at the incompressible insulator-like 
Cu sites. Inspection of the sum rule in equation (|5| shows that the incompressible Cu3<i 
states and the compressible 02p states are entangled by the off-diagonal elements of 
the polarizability matrix. So our modelling of the pseudogap state is consistent with a 
strange not overall compressible metallic state coexisting with atom-like, incompressible 
Cu3d orbitals similar as in the undoped, insulating state favoring short ranged magnetic 
spin correlations and corresponding excitations, typically spin fluctuations. The latter 
should be damped by electron-hole charge excitations of the compressible metallic 02p 
states. 

In /c-space language such a strange mixed locally compressible-incompressible state 
should result in a heavily renormalized electronic bandstructure and Fermi surface (FS) 
as compared with the more conventional overall compressible Fermi liquid state with a 
large FS at higher doping. The pseudogap state displaying properties characteristic of 
both, insulators and correlated metals evolves from the incompressible insulating state 
with no low energy charge excitations and seems to promote the tendency to charge 
and spin ordering that breaks rotational and translational symmetry, e.g. in form of 
stripes. The latter appear to be quite common in underdoped cuprates. In this context 
it is interesting to point out that the stripe-like modulation observed in [5] is strongest 
on the oxygen orbitals. So, the latter are very essential to understand the physics of 
charge order in the cuprates and also the charge response as shown below. Finally, 
the enhanced phase fluctuations due to the incompressibility of the Cn3d orbital favour 
pairing correlations in the normal state, possibly incoherent preformed pairs as discussed 
in the literature. 
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(a)Experimental Results (b)Calculated Results 

Figure 1. (a) Experimental results for the highest Ai and Si branch of La2- a; Sr 2 ;Cu04 
for various doping levels (Ref. 4). •: x = 0.00, ■ : x = 0.10, o: x = 0.15, □: 
x = 0.30. The lines are a guide to the eye. (b) Calculated results of the phonon 
branches shown in (a) as obtained with our model approach for the electronic state of 

the HTSC's [T]. : insulating state, underdoped state, : optimally doped 

state, : overdoped state. For comparison, the results for the RIM are shown ( — ) 

in order to demonstrate the crucial influence of the nonlocal EPI effects in form of CF's 
on the dispersion. 

Our modeling of the density response is strongly supported by corresponding 
calculations of phonon dynamics [JJ, [2] , which compare well with the experimental data. 
This is shown for p-doped LaCuO in case of the strongly doping dependent OBSM 
phonon anomalies, see figure [TJ where the experimental results for the highest Ax and 
Sx branches of L&2- x Sy x CvlO<i are shown for various doping levels [6] and compared 
with the calculated results. The high-frequency OBSM are generally known to display an 
anomalous softening upon doping in the cuprates [ZllHlElIinillllinillSlElIISlIISlIIZllIH]. 
Such a behaviour of the experimental phonon dispersion is clearly visible for LaCuO in 
figure [T^i and is very well accounted for by the calculation in our model for the electronic 
state in terms of incompressibility-compressibility transitions, see figure [T]d. We also 
have included in figure [T]d the phonon dispersion of a rigid-ion model (RIM) for LaCuO. 

It seems appropriate to explain at least the ideas of our modeling of phonon 
dynamics and charge response of the cuprates. A quantitative treatment of the theory 
and the modeling would be too extensive and can be found in Refs. [TJ [2]. 
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The rigid part of the electronic charge response and the EPI is approximated 
by an ab initio rigid-ion model (RIM), taking into account covalent ion softening in 
terms of (static) effective ionic charges calculated from a tight-binding analysis. In 
addition, scaling of the short range part of certain pair potentials between the ions 
is performed to simulate further covalence effects in the calculation in such a way 
that the energy-minimized structure is as close as possible to the experimental one. 
The RIM with the corrections just mentioned then serves as an unbiased reference 
system for the description of the cuprate superconductors and can be considered as a 
first approximation for the insulating state of these compounds because of the strong 
ionic nature of bonding in the cuprates. Starting with such a rigid reference system 
nonlocal, electronic polarisation processes are introduced in the form of more or less 
localized electronic charge fluctuations (CF's) at the outer shells of the ions. Especially 
in the metallic state of the cuprates the latter dominate the nonlocal contribution of 
the charge response and the EPI and are particularly important in the CuO planes. In 
addition, anisotropic dipole fluctuations are admitted in our approach. They prove to 
be specifically of interest for the ions in the ionic layers. 

Comparing in figure [T]d the results of the phonon dispersion of the RIM, where 
only local EPI effects are considered, with the full calculation we can extract the crucial 
influence of the nonlocal EPI effects generated by the charge fluctuations. For example 
we find a very strong softening of the Ai/2 bond-stretching mode (half-breathing mode) 
increasing with doping up to about 5 THz. Comparing the calculated results with 
the experiments in figure [T^i we find a good agreement. In order to achieve such a 
realistic description of the phonon modes, both, nonlocal EPI effects and multi-orbital 
effects are required. The comparison with the results from the RIM demonstrates that 
nonlocality is needed and from our detailed calculations we find that the CF's of the 
more delocalized 02p and Cu4s orbital contribute significantly to the softening of the 
modes, besides the more localized Cu3<i orbital. The strong softening of the phonon 
anomalies is accompanied by corresponding dynamic charge inhomogenities in form of 
charge stripes along the CuO bonds. The modes are not really soft but one may consider 
these dynamic stripes as precursors of static charge stripe order with ordering vectors 
corresponding to the wavevectors of the phonon anomalies. 

Finally, it should be remarked that a constructive interplay of electron-ion 
interaction with short-ranged antiferromagnetic order is important for the Cu sites in 
the underdoped and insulating state, because electrons are transferred from the Cu ion 
where the Cu bond is compressed to the Cu ion where the bond is stretched, see e.g. 
figure [2j Both ions have opposite spin in the antiferromagnetic case. In contrast if a 
ferromagnetic order would be present the Pauli blocking would suppress the transport 
of electrons with the same spin and consequently also the charge fluctuations and the 
corresponding softening of the phonon modes. 

Modelling of the electronic polarizability H KK i (q ) for the insulating and pseudogap 
state is needed because a quantitative calculation of IT reK / (q ) is not available due 
to the missing of a rigorous knowledge of this state. On the other hand, for 
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Figure 2. (a) Contour plot in the CuO plane of the displacement induced charge 
density redistribution 5p for the Ai/2 mode within the rigid ion model (RIM) for 
La2CuC>4. (b) Same as (a), but for the nonlocal part of the charge response in the 
optimally doped state due to nonlocal EPI effects in the half-breathing mode Ai/2. 

The displaced ions are indicated by a dot with an arrow. Full lines ( ) mean that 

electrons are accumulated in the corresponding region of space and the broken dotted 

lines ( ) indicate regions where the electrons are pushed away. Sp is given in units 

of 10- 4 e 2 /a|- 



the more conventional metallic state H KK ' (q ) can be calculated within the local 
density approximation (LDA), however, with some important corrections concerning 
the electronic c-axis response which is overestimated in typical LDA-based calculations 
[19J. Concerning this case it should be remarked that a comparison of the calculated 
OBSM anomalies from the model in figure [T]d with corresponding results obtained within 
the framework of the corrected LDA [19], taking into account the full q dependence of 
the polarizability matrix, leads to a good agreement. This also speaks in favour of our 
modeling of the charge response in terms of the compressibility sum rules, in which only 
the longwavelength limit of the polarizability is considered. 

In figure [2] we display the calculated phonon induced charge redistribution of the 
Ai/2 bond-stretching mode for the optimal doped state. Figure [2^, illustrates the 
effect of the local EPI within the RIM connected with the displacement of the O x 
ions in the CuO plane, i.e. a stretching or compression of the CuO bond, respectively, 
(half-breathing mode (HBM)). The local EPI reveals a charge redistribution of dipole 
character localized only at the ions as can be expected from the rigid shift of the ionic 
densities in the RIM during the phonon mode. 

From figure [2^, we also obtain the message that no dynamical charge stripes are 
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(a) (b) 

Figure 3. (a) Same as figure 2 but for the nonlocal part of the charge response Sp for 
Ai/2 in the insulating state (x = 0). (b) Same as figure 2 but for the nonlocal part of 
the charge response 5p for Ai/2 in the pseudogap state (x = 0.1). 



excited by local EPI. On the other hand, figure^ that displays the nonlocal part of the 
charge response due to nonlocal EPI effects of charge fluctuation type demonstrates, 
how the strong nonlocal EPI generates dynamic charge stripes in the half-breathing 
mode. Thus, the moving O x ions in the HBM which are linked up with strong CF's on 
the silent Cu and O y ions induce dynamical stripes in the y direction with a recurrence 
periond of two lattice constants along the orthogonal x-direction. Note, that there are 
no changes of the transfer integral between d and p orbitals for the silent Cu and O y 
ions, nevertheless there is a charge transfer. 

In case of the full breathing mode Og, at the the X point of the BZ where the O x 
and O y ions vibrate in phase along the bonds against the silent Cu ion pQ nonlocal EPI 
generates CF's at the Cu and we obtain an electronic charge transfer from that Cu ion 
where the bonds are compressed to the Cu where the bonds are stretched [20J. While 
in case of the HBM the charge stripes point along the x- or ^/-direction, respectively, for 

the stripes point along the diagonals in the CuO plane. 

The self-consistent changes of the potential felt by an electron due to nonlocal 
EPI also have been calculated for the OBSM and yield a strong coupling of the order 
100 meV [IH]. It should be remarked that such a strong nonlocal EPI associated with 
the OBSM, which can reliably be calculated in the adiabatic approximation, is further 
enhanced by an order of magnitude because of the poorly screened long range Coulomb 
interaction and the slow charge dynamics for the polar modes from the nonadiabatic 
sector of phase space around the c-axis where phonon-plasmon mixing is relevant |19j . 
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In figure |3^ and figure [3]d we present the calculated results for the nonlocally 
induced charge redistribution of the HBM for the insulating and the pseudogap state, 
which should be compared to the corresponding results for the optimal doped state 
in figure [2Jd. First of all we find for all electronic states dynamical stripes along the 
bond direction excited by the nonlocal EPI being simultaneously responsible for the 
corresponding phonon softening in figure [TJ This is quite interesting because dynamic 
charge inhomogenities of this type may be considered as precursors of static charge 
stripe order as recently observed in L&2- x Ba x . CuO/i, in a broad range of doping around 
x = 1/8 with 0.095 < x < 0.155 [21J or in YBa 2 Cu30 y , respectively, [22]. 

To the best of our knowledge in overdoped samples static stripe order has not 
been observed so far. This would mean that the pseudogap in the underdoped state 
may be required to nucleate dynamic stripes. In our approach this implies that for 
the emergence of static stripe order the incompressibility of the Cu3d orbitals seems 
to be essential. The latter brings about a separation of mobile holes populating the 
compressible 02p states and regions with local Mott- and spin correlations related to 
the insulator like Cu3d states. Such stripes producing the onset of a new periodicity that 
breaks translational and rotational symmetry also may help to trigger a reconstruction of 
the large hole Fermi surface (FS) into small pockets as detected by quantum oscillation 
measurements in underdoped YBa 2 Cu 3 6+x , see e.g. [23] and references therein. Also 
from an energetic point of view our modeling of the pseudogap state is susceptible for 
the tendency to establish broken symmetry states with hole-rich compressible and hole- 
poor incompressible regions. In general terms lowering of the kinetic energy occurs in 
the compressible metallic regions and lowering of the interaction energy is possible by 
antiferromagnetic spin correlations and the Coulomb repulsion between the electrons in 
the region of the incompressible insulator-like Cu3d states. Altogether, in the pseudogap 
state the doped holes and the superexchange-coupled spins arrange themselves in 
a cooperative way to promote both charge mobility and locally antiferromagnetic 
correlations. 

Comparing the stripe patterns of the insulating and underdoped state with that 
of the more conventional metallic optimal and overdoped state (not shown) with a 
large FS as the locus of the crucial low energy excitations, we find a relatively small 
contribution of the 02p orbitals degrees of freedom to the density redistribution in the 
latter case. The weight of the 02p contribution is significantly enhanced in the insulating 
and pseudogap state, see figure|3^ and figure^ reflecting in particular the importance of 
the 02p orbital in the pseudogap state. Holes in the pseudogap state primarily populate 
the 02p orbitals, while the holes in the more conventional metallic state gain more and 
more Cu3d character by an increased hybridization as can be seen from a large PDOS 
for Cu3<i leading to a dominating contribution of the charge redistribution at the Cu 
ion where the charge response is now metallic, too. See also [24] where it is investigated 
how the nature of the hole state evolves with doping. Overall, the stripe pattern is more 
delocalized in the metallic state and localization increases towards the insulating state 
via the pseudogap state. 
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Figure 4. Same as figure 2 but for the nonlocal part of the charge response Sp for 
Ai/4 in the optimally doped state. 

In the overdoped state a growing importance of the Cu4s orbital and accordingly an 
increase of the Cu4s element of the polarizability matrix leads to an enhanced softening 
of the anomalous OBSM, see figure [T] The corresponding calculated dynamical stripe 
pattern is very similar to that of the optimally doped state in figure [2]d but somewhat 
more delocalized because of the increased contribution of the more delocalized Cu4s 
contribution. 

From our calculations it becomes obvious that a multi-orbital approach is needed 
that includes besides Cn3d at least 02p orbitals and for the optimally doped and 
overdoped state in addition Cu4s orbitals in order to give a reliable representation of 
phonon dynamics and charge response in the cuprates. Single-band models for the 
CuO plane believed by many people to describe the cuprates are not sufficient, in 
particular in context with the discussion of stripes. Altogether, in the optimal and 
overdoped state increasing of doping enhances the hybridization between Cu3<i, Cu4s 
and 02p orbitals, respectively, with a dominant weight for Cu3gL This state becomes 
compressible metallic with a large PDOS at ep. So, the particle number fluctuations 
that characterize the different competing ground states in our model are enhanced as 
compared to the ground state in the pseudogap regime. Here the hybridization is weaker 
due to the incompressible insulator-like Cu3d orbital with a vanishing PDOS at ep. 

Translational symmetry breaking models of the electronic structure of the CuO 
plane generated by charge density wave order of a certain wavevector are discussed 
in the literature and used to characterize Fermi surface reconstructions of a large FS 
into small pockets, for recent work see f2~5[ [2S]- From scanning tunneling microscopy 
(STM), neutron diffraction experiments and nuclear magnetic resonance measurements 
the ordering in many cases occurs with a period of about 4 lattice constants [26]. This 
corresponds to the wavevector j- (1, 0, 0) of the Ax/4 mode. Therefore it is meaningful 
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to calculate the nonlocally induced charge redistribution also for this mode. The result 
for the optimal doped state is shown in figure |4| 

We again obtain a stripe like charge pattern with translational symmetry in the 
(0,1,0) direction but a four unit cell repeat distance along the orthogonal direction 
(1,0,0). Contrasting this result with the charge stripe pattern as found for the HBM in 
figure [2}} the stripes for Ai/4 are expanded over the whole unit cell, because the CF's do 
not vanish at the sites of the moving oxygens, unlike to Ax/2. However, the amplitudes 
of the CF's are considerably smaller. According to our calculation of the eigenvectors 
this is caused by the strongly reduced bond-stretching amplitude of the oxygen ions in 
the plane for Ai/4, while the displacement of the apex oxygen ions is strongly enhanced 
in the Ai/4 mode as compared with the HBM. 

In summary, we have applied our modeling of the doping dependent electronic 
states of the cuprates to the calculation of the phonon induced charge inhomogeneities 
nonlocally induced by the Ax/2 and Ax/4 mode in LaCuO. The calculations have been 
performed for the insulating-, pseudogap and the more conventional metallic state. The 
compressibility is used as a primary term to characterize the competing ground states 
controlled by doping. The strong nonlocal EPI of CF type excites in all cases investigated 
dynamical charge stripes along the CuO bonds. Differences in orbital character of 
the charge response in the various states are discussed at hand of the stripe patterns. 
Experimentally, charge stripe order has been observed in a broad range of doping in 
La2_ a; Ba 3 ;Cu04 an d other cuprates which supports the point of view that dynamical 
charge stripes generated by certain strongly coupling OBSM may be considered as 
precursors of static stripe order at least outside the overdoped regime. Such a charge 
order breaks translational symmetry and may enforce a reconstruction of a larger hole 
FS into small pockets. 

The interplay between the reduced particle density fluctuations related to 
the insulator-like behaviour of the localized atomic-like Cu3d orbital that remains 
incompressible for low enough p-doping and on the other hand the compressible, metallic 
02p orbitals appear to play an important role for the pseudogap state and the tendency 
to charge ordering in the cuprates. 
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